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We study the resonant x-ray scattering (RXS) spectra at the Dy Lni absorption edge in the 
quadrupole ordering phase of DyB2C2. Analyzing the buckhng of sheets of B and C atoms, 
we construct an effective model that the crystal field is acting on the 5d and 4/ states with 
the principal axes different for different sublattices. Treating the 5d states as a band and the 
4/ states as localized states, we calculate the spectra within the dipole transition. We take 
account of processes that (1) the lattice distortion directly modulates the 5d states and (2) the 
charge anisotropy of the quadrupole ordering 4/ states modulates the 5d states through the 
5d-4/ Coulomb interaction. Both processes give rise to the RXS intensities on (00 1) and (/iO|) 
spots. Both give similar photon-energy dependences and the same azimuthal-angle dependences 
for the main peak, in agreement with the experiment. The first process is found to give the 
intensities much larger than the second one in a wide parameter range of crystal field. This 
suggests that the main-peak of the RXS spectra is not a direct reflection of the quadrupole 
order but mainly controlled by the lattice distortion. 
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1. Introduction 

\ Resonant x-ray scattering has recently attracted much 

■ interest, since the resonant enhancement for the pro- 
\ hibited Bragg reflection corresponding to the orbital or- 

■ der has been observed in several transition-metal com- 
\ pounds by using synchrotron radiation with photon en- 

■ ergy around the K absorption edge.^"''^ For such K- 
\ edge resonances, 4p states of transition metals are in- 

■ volved in the intermediate state in the electric dipolar 
\ [El) process, and they have to be modulated in accor- 

■ dance with the orbital order for the signal to be observed. 
\ This modulation was first considered to come from the 

■ anisotropic term of the 4p-3(i intra-atomic Coulomb in- 
\ teraction,^^ but subsequent studies based on the band 

■ structure calculation^"^^ have revealed that the modu- 
I lation comes mainly from the crystal distortion via the 

■ oxygen potential on the neighboring sites. This is be- 
\ cause Ap states are so extending in space that they are 
i very sensitive to the electronic structure at neighboring 
' sites. 

Rare-earth compounds also show the orbital order 
(usually an ordering of quadrupole moments). In CeBg, 

■ RXS experiments were carried out around the Ce Lm 
' absorption edge, and resonant enhancements have been 

found on quadrupolar ordering superlattice spots. 
Only one peak appeared as a function of photon energy, 
which was assigned to the Ei process. In the Ei process, 
5d states of Ce in the intermediate state are to be mod- 
ulated in accordance with the superlattice spots. Since 
the lattice distortion seems extremely small and the Sd 
states are less extending than the Ap states in transition- 
metal compounds, it is highly possible that the modu- 
lation is mainly caused by the Coulomb interaction be- 
tween the bd states and the orbital ordering 4/ states. 
In our previous papers, we demonstrated this sce- 
nario by calculating the RXS spectra on the basis of the 
effective Hamiltonian of Shiina et al.^'^"^'^-' Without the 



help of lattice distortion, we obtained sufficient intensi- 
ties of the spectra, and reproduced well the temperature 
and magnetic field dependences. This situation contrasts 
with those in transition-metal compounds. 

Another example for rare-earth compounds is RXS ex- 
periments on DyB2C2, where the intensity is resonantly 
enhanced near the Dy Lm absorption edge.^^"^^^ This 
material takes a tetragonal form at high temperatures 
as shown in Fig. 1(a), and undergoes two phase transi- 
tions with decreasing temperatures in the absence of the 
magnetic field: a quadrupole order below Tq (= 24.7 K) 
(Phase II) and a magnetic order below Tq {— 15.3 K) 
(Phase III).^^) Corresponding to the transition at Tq, 
a large non-resonant intensity is found in the a a' 
channel on the (/iO|) spot {h and i are odd integers). 
This suggests that some structural change takes place 
at T = Tq from the tetragonal phase at high temper- 
atures.^^' A buckling of sheets of B and C atoms 
was proposed, and the non-resonant intensities by the 
buckling has recently been evaluated; about 0.01 A shift 
of B and/or C atoms may be sufficient to give rise to 
such large intensities.^*'-' It is not clear in experiments 
whether the intensity on this spot is resonantly enhanced 
at the Lm edge, since the non-resonant part is so large 
that it may mask the resonant behavior. On the other 
hand, the resonant enhancement of RXS intensities has 
clearly been observed on the superlattice spot (00|). 

In this paper, we study the mechanism of the RXS 
spectra at the Lm edge in Phase II of DyB2C2. Since 
the 5(i states are so extended in space that they are sensi- 
tive to lattice distortion caused by the buckling of sheets 
of B and C atoms. Then the question arises whether the 
direct influence of the lattice distortion on the 5(i states is 
larger than the influence of the anisotropic 4/ charge dis- 
tribution associated with the quadrupole order through 
the 5d-4/ Coulomb interaction. Lovesey and Knight^^^ 
have discussed the mechanism from the symmetry view- 
point, and have pointed out that the RXS intensities on 
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Fig. 1. (a) Sketch of the crystal structure of DyB2C2 {PA/mbm: 
a = 5.341 A, c = 3.547 A at 30 K). Gray large circles are Dy 
atoms. Solid and open small circles are B and C atoms, respec- 
tively, (b) Local coordinate frames attached to each sublattice. 



Note that the Ce-Ce interaction in CeBg has been exten- 
sively studied, describing well the phase diagram under 
the magnetic field. But the molecular field may 
change little and even stabilize the quadrupole order. 
Therefore, we need not explicitly consider the molecular 
field by regarding the crystal field as including the effect. 
The charge anisotropy associated with the quadrupole 
order modulates the 5d states through the intra-atomic 
5d-Af Coulomb interaction. 

We calculate the RXS intensity within the Ei transi- 
tion. We take account of the above two processes, direct 
and indirect ones, of modulating the 5d states. Both 
processes give rise to the RXS intensities on the (00 1) 
and on the (ft.O|) spots. Both give similar photon-energy 
dependences and the same azimuthal-angle dependence 
in agreement with the experiment. However, the mech- 
anism of direct modulation of the 5d band gives rise to 
the intensities much larger than the mechanism of indi- 
rect modulation through the 5o?-4/ Coulomb interaction 
in a wide parameter range of the crystal field. This sug- 
gests that the RXS intensities are mainly controlled by 
the lattice distortion. 

This paper is organized as follows. In § 2, we analyze 
the buckling of sheets of B and C atoms. In § 3, we 
briefly summarize the formulae used in the calculation of 
the RXS spectra. In § 4, we calculate the RXS spectra 
on two mechanisms. Section 5 is devoted to concluding 
remarks. 



(00|) and (ft-0|) spots come from lowering the local sym- 
metry probably due to lattice distortion. The argument 
based on symmetry alone is powerful in some respect, but 
does not shed light on this issue. In the transition-metal 
compounds, the corresponding question has already been 
answered by ah initio calculations as mentioned above. 
However, such ah initio calculations are difficult in rare- 
earth compounds. We resort to a model calculation by 
treating the 5d states as a band and the 4/ states as lo- 
calized states. The buckling of sheets of B and C atoms 
causes modulations of the 5c? bands and of the 4/ states. 
We analyze such effects of lattice distortion on the ba- 
sis of the point charge model, which leads to four in- 
equivalent Dy sites with principal axes of the crystal field 
shown in Fig. 1(b). These principal axes seem to corre- 
spond well to the direction of magnetic moments in the 
magnetic phase. Of course, the point charge model 
is not good in quantitative viewpoint. Nonetheless, we 
construct an effective model that the 5d and 4/ states 
are under the crystal field of the same form and with the 
same principal axes as the above analysis. 

The crystal field modulates the 5d states. Although 
the actual effect may come from hybridizations to 2p, 
3s states of B and C, it can be included into a form 
of the crystal field. The crystal field also makes the 
quadrupole moment of the 4/ states align along the prin- 
cipal axes, establishing a quadrupole order. A molecular 
field caused by the Dy-Dy interaction may also act on 
the 4/ states in Phase II in addition to the crystal field. 
This interaction may be mediated by the RKKY inter- 
action, but the explicit form has not been derived yet. 



2. Lattice Distortion 

For making clear the effect of lattice distortion on elec- 
tronic states, we first calculate the electrostatic potential 
on the basis of a point charge model. Point charges 
iDy, q_B, and qc are placed on Dy, B, and C sites, respec- 
tively. Figure 2 shows the positions of B and C atoms on 
the plane of z = c/2 and those of Dy atoms on the plane 
of z = 0. For the buckling of sheets of B and C atoms 
(up- and down-movements along the c axis) specified in 
Fig. 2, the electrostatic potential is evaluated within the 
second order of coordinates around Dy sites. As shown 
in Fig. 1(b), four inequivalent sites arises. For sites j 
(— 1 • • • 4), the electrostatic potential is obtained as 



VcrysU) = AlQl + Al{j)Ql + A,y(j)Q^y, (1) 



with 



Qxy = V^xy, 
where the coefficients are given by 
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z = 1/2 




6lg) = 180° - 65.6°, e^c^ = 180° - 19.1°, 

= 90° + e^g^ , 6»^^^ = 90° + e^c^ , 



90° + 6'[?\ 



(10) 



Now we search for the local coordinate frames in which 
the third term in eq. (1) is eliminated. Rotating the 

original coordinate frame by angle around the c axis 
for each sublattice, we have the operators transformed 
as 

Ql = cos{2(pj)Ql{j) - sm{2(f>j)Q^y{j), 

Q„j = sin(2<^,)Q2(j) + cos(20,)Qxv(i), (H) 

where tilde operators Q(j)'s are represented with respect 
to the local coordinate frames. Q2 is unchanged. Insert- 
ing eq. (11) into eq. (1), we have 

Vorysij) = AlQ%{j) + Al{j)Ql{j) + A^y{j)Q^y{j), (12) 




Dy 



Fig. 2. Sketch of a B2C2 sheet {z = c/2). Open circles represent 
B and C atoms; big and small circles move to positive and nega- 
tive directions along the z axis, respectively. The directions are 
reversed on the plane oi z = —c/2. Hatched circles represent Dy 
atoms at the plane of 2 = 0. 



2qD 



l + 1.4oM + i.79l^c;| 



(ID 



QD 



Axyij) 
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The first term in eq. (1) represents the crystal field with- 
out lattice distortion, while the second and third terms 
arise from the buckling. The Rb (= 2.732A) and Rc 
(= 2.676A) are distances from the origin to B and C 
sites, respectively. The dB and dc represent the abso- 
lute values of shifts along the c axis from the z = ±c/2 



planes, respectively. Angles ^^(C) ^s(C) (^) 
and (7) are given by 

3(1) 



0(1) 



65.6°, e 



19.1° 



with 



Ai{j) = Ab cos[2(eg) - cPj)] + Accos[2(e^^^ - </,,)], 



(13) 



ixy(j) = ABsin[2(0(^^ 



Acsin[2(0g^ -0,)] 



(14) 



Condition A^yij) = determines (pj^s, which take val- 
ues between 0^^ and O'^K For example, assuming qb = 
qc = — (3/4)e, qoy = 3e (c: proton charge), dB = O.OlA, 
dc = 0.02 A, we have = 31.8°, 02 = 180° - (j)!, 
03 = 90° -I- 01, (/)4 = 90° -|- 02- The principal axes thus 
estimated correspond well with the directions of the or- 
dered magnetic moments in Phase III.^^^ 

The equivalent operator method allows us to write the 
crystal field energy Hcrys (j) at site j within the subspace 
of angular momentum J as 

H^rysU) = DjWl{j) - J{J + 1)] + Ej[jl{j) - Jim. 

(15) 

(a) on the 5d bands 

The 5d states are forming an energy band with width 
^ 15 cV through a hybridization with s and p states 
of neighboring B and C atoms as well as 5d states of 
neighboring Dy atoms. We need the density of states 
(DOS) for calculating the RXS intensity. We assume a 
Lorentzian with full width of half max;imum 5 eV for 
the DOS's projected onto symmetries xy, — y^, y^;. 



2X, and Jtz^ 



The center of each component of the 



DOS is separate to each other in acc;ordance with the 
first term of eq. (15), although the first term need not 
be respected so much because of the large band efii'ect. 
Explicitly they are assumed to be (l/7r)A/((e — 2.5)^ -|- 



A^) for xy. 



y\ (lA)A/((e-7)2 + A2) for yz, zx, 



and (l/7r)A/((e-8.5)2-^A2) for S^^-r^, with energies in 
units of cV and A = 2.5 eV. This arbitrary assumption 
for the DOS form may be justified by the fact that the 
RXS spectra is not sensitive to the assumption. 

The second term of cq. (15), which arises from the 
buckling of sheets of B and C atoms, gives rise to a 
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small modification on the 5d band. Although the ac- 
tual modulation of the 5d states may come through the 
hybridization to the s and p states of B and C atoms, 
such effects can be includcid into the second term. This 
term makes the local symmetry twofold. 
(b) on the 4f states 

Dy"^"*" ion is approximately in the 4/^-configuration 
(^Hi5/2). Equation (15) is now applied to the subspace 
of J = 15/2. Since the 4/ states are much localized than 
the 5(i states, the crystal field is much smaller here than 
that on the 5d states. The coefficient -D/ of the first 
term is expected to be positive from the analysis of mag- 
netic susceptibility.^^) This leads to the lowest energy 
states I ± ^) and the next energy states | ± |), both of 
which form Kramers' doublets (|M) represents the state 
of Jz = M). The axial symmetry is kept instead of four- 
fold symmetry without the lattice distortion. Terms of 
of = + J1), which admixes the states \M) with 

|M ± 4), come from the higher order expansion to make 
the local symmetry fourfold. The detailed study along 
this line, however, is beyond the scope of the present 
study. 

In any event, the second term makes the local sym- 
metry twofold. The lowest energy state is admixed by 
|M) with \M\ > 1/2. The quadrupole moment is or- 
dered; {0,j.2_y2) (= ^{Jx^Jy}) takes a finite value with 

in the local coordinate 
) means the average over 
the lowest Kramers doublet). With increasing values of 
\Ef\/Df, the average quadrupole moment increases, as 
shown in Fig. 6(b). It becomes largest, {0^-2_y2) = 46.1 
at Df = 0.1. 

3. Cross Section of Resoneint X-Ray Scattering 

We briefly summarize here the formulae used in the 
calculation of the RXS spcx'tra in the next section. The 
conventional RXS geometry is shown in Fig. 3; photon 
with frequency w, momentum kj and polarization ^J. {= a 
or it) is scattered into the state with momentum ky and 
polarization fj,' (= a' or tt'). The scattering vector is 
defined as G = k/ — kj. Near the Dy Lm absorption 
ed ge, a 2p core electron is virtually excited to 5d states 
in the Ei process. Subsequently it recombines with the 
core hole. Since the 2p states are well localized around 
Dy sites, it is a good approximation to describe the scat- 
tering tensor as a sum of contributions from each site 
of the core hole. Therefore, the cross section in the Ei 
process is given by 



frame for each sublattice ((• ■ 



(Oxy) (— ^2 ^JxJy 



(16) 



where 



.(G,a;) = --^^M„„,(i,w)exp(-iG-r,), (17) 



with 

Maa' {j, W) 



E 



(j)|A)(A|x^,(j)|^,-) 
hio - {Ea - Ej) + iT '■ 



(18) 



Table I. Geometrical factors 



1 COS /3 cos tp cos /9 cos "0 

2 — sin i/i — sin i/> 

3 — sin /3 cos V' — sin /3 cos V' 



- sin 9 cos (3 sin ip + cos sin /? 
— sin 9 cos i/' -l- cos 9 sin jS 
cos 6 cos /3 




Fig. 3. Scattering geometry. Incident photon with wave vector 
kj and polarization it or tt is scattered into the state with wave 
vector kj and polarization cr' or tt' at Bragg angle 9. The sample 
crystal is rotated by azimuthal angle V around the scattering 
vector G = kj: — k,. 



where N is the number of Dy sites. Note that the 
cross section is an order of A^. The P'' and P'^^ are 
geometrical factors for the incident and scattered pho- 
tons, respectively. Their explicit forms are given in Ta- 
ble I. The represents the initial state with energy 
Ej. The intermediate state |A) consists of an excited 
electron on 5c? states and a hole on 2p states with en- 
ergy Ef^. The F is the life-time broadening width of 
the core hole. The dipole operators XaijYs are defined 
as xi{j) = .T, X2{j) = y, and x^{j) = 2; in the coor- 
dinate frame fixed to the crystal axes with the origin 
located at the center of site j. The scattering ampli- 
tude Maa'{G,uj) contains the square of the dipole ma- 
trix clement Adp — (5d|r|2p) = R5d{r)rR2p{r)r'^dr 
with R5d{r) and R2p{r) being the radial wavefunctions 
for the 5d and 2p states. For Dy^+ atom, it is estimated 
as 2.97 X 10~^^ cm in the 4/^ configuration within the 
HF approximation.^'^) 



4. Calculation of RXS spectra 

In the Ex transition, an electron is excited from 2p 
states to 5d states at a Dy site. A 2p core hole is left 
behind, and its state is split into the states of jp = 3/2 
and jp = 1/2 (jp is the total angular momentum) due 
to the strong spin-orbit interaction. We consider only 
the jp = 3/2 states (Lm edge). We describe the pho- 
toexcited bd electron in the band by introducing a local 
Green's function. 



g: 



5d 



iS 



de. 



(19) 



where p^<i(e) is the m'^ component of the d DOS defined 
in §2. The Fermi energy e_F is set to be zero so that the 
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Table II. Slater integrals and the spin-orbit intera<;tion for Dy^+ 
atoms in the Hartree-Fock approximation (in units of eV). 



F>'(4f,4f) 
F° 32.19 
15.31 
9.607 
F6 6.911 



F>'{2p, 5d) 
F° 16.13 
0.489 



F*(2p,4/) 
F° 44.62 
1.982 



F° 14.70 
F2 3.614 
1.741 



G*(2p,5d) 
Gi 0.414 
G3 0.245 



G*(2p,4/) G'=(4/,5(i) 
G2 0.207 Gi 1.615 
G^ 0.133 G^ 1.321 
G^ 1.009 



C4/ = 0.273 Csd = 0.181 



*In the RXS calculation, the above values of the anisotropic terms 
axe reduced by multiplying a factor 0.8, while the values for 
F^^\nl,n'l') are replaced by much smaller values, F^") (4/, Sd) = 
3.0, F(o)(4/,4/) = 7.0, F(0)(2p,5d) = 4.0, F(o)(2p,4/) = 12.0. 



5d band is almost empty. 

Finite RXS intensities on superlattice spots arise from 
modulating the 5d states with wave vectors of superlat- 
tice spots. There are two origins to giving rise to such 
modulation. One is a direct influence from the buckling 
of sheets of B and C atoms, which is represented by the 
second term of eq. (15). Another is the charge anisotropy 
of the 4/ states in the quadrupole ordering phase. We 
discuss both origins separately. In the actual calculation, 
we specify the local coordinate frames with (jji = 28°, 
(/)2 = 180° - 28°, (/)3 = 90° + 28°, (/)4 = 90° + 152° for 
four sublattices in accordance with the experiment. 

4-1 Direct influence of lattice distortion 

Let the Ei transition take place at a particular site 
(called as "origin"). The excited 5d electron is attracted 
by the core hole potential at the origin. What is more 
important is that the 5d electron is under the influence 
of the second term of eq. (15). Taking account of the 
multiple scattering from these terms at the origin, we 
evaluate the resolvent l/{hu! — Hint) with respect to the 
intermediate-state Hamiltonian Hi^t- 

1 



hu) — Hu 



iT 



(20) 



where m'' and s'^ specify the orbital and spin of the d 

electron, respectively. The ex represents the energy of 
the core hole with A in the jp = 3/2 subspace. The 
scattering potential Vfy^dsdx-.m'ds'dX' includes the second 
term of eq. (15) and the Coulomb interaction between 
the 5d electron and the 2p hole. The latter quantity 
is expressed in terms of the Slater integrals, which are 
evaluated within the HF approximation in a Dy'^+ atom 
(see Table II). The core hole life-time width is set to 
be r = 2.5 eV. Equation (20) is numerically evaluated. 

Before calculating the RXS spectra, we touch on the 
absorption coefficient. We calculate the absorption co- 
efiicient A{w) in the Ei process from the resolvent by 
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Exp. 
T=7K 
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Absorption Coeff. 



7790 



7810 



a)[eV] 



Fig. 4. Absorption coefficient A{ui) (lower panel) in comparison 
with the Ljii-edge fluorescence spectra (upper panel). 



using the relation. 



A(u;)(x^^(^,>"(j)|mVA) 

j a 



X I — — ) Im ( 

TtJ - Hint + iO / ^dsdX;m"is"iX' 



X (m"*s«A>«(i)|^,-). 



(21) 



Here ImX indicates the imaginary part of the quantity 
X. Figure 4 shows the calculated A{u)) in comparison 
with the fluorescence experiment. We adjust the core 
hole energy such that the calculated peak coincides with 
the experimental one. The calculated curve reproduces 
the experimental one. Since the spectra is proportional 
to the d DOS when the 5d-2p Coulomb interaction is 
neglected, the assumed DOS seems reasonable. The in- 
tensity seen in the high energy region of the experiment 
may come from the d symmetric states mixing with 3s, 
3p states of B and C atoms, which is outside our interest. 

Now we discuss the RXS spectra. The resolvent above 
calculated is used to calculate the scattering amplitude 
at the origin with the help of the relation, 



E 



|A)(A| 



huj - {Ea - Ej) + iT 



E E 



1 



hjjo — Hi, 



■iT 



{m"'s"^X'\. (22) 



This expression is independent of the quadrupole order- 
ing 4/ states. It is inserted into eq. (18) to calculate 
an RXS amplitude. The extension to general site j is 
straightforward. In the coordinate frame fixed to crys- 
tal (not in the local coordinate frames), their forms are 
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7(0;) 



C(w) 

-6(uj) 

7(0;) 

^(w) -^(w) 

C(w) 

7(w) 



(23) 



(a) G = (00|) 7OTf/i ^ 0(i(i integers. 

Owing to the factor exp(— zG • r^) in eq. (17), the 
total amplitude is given by a combination of M(l,w) + 
M(2,cj) - M(3,a;) - M(4,a;). Thus, we have the final 
form, 



M(G,w) 



N 




(24) 



The geometrical factors are given by setting [3 — in 
Table I, which are combined to eq. (24) to calculate 
the scattering intensity. We have the RXS intensity as a 
function of azimuthal angle ip as 

oc |77(a;)|^sin^2V', 

7^^„/(G,a;)oc |77(a;)|^ sin^ ^cos^ 2^-, (25) 

with 9 the Bragg angle. Here ip is defined such that 
tp = corresponds to the scattering plane containing the 
b axis. 

Figure 5 shows the calculated RXS spectra as a func- 
tion of photon energy in comparison with the experiment 
(t/) = —45°).^^^ The crystal field parameter is set to be 
Ed = —0.1 eV in eq. (15). As shown in the middle panel, 
the calculated spectra show a single-peak in agreement 
with the experiment. (Only the a ^ a' channel gives 
finite intensity for i/j = —45°.) The photon energy de- 
pendence in the ct — > tt' channel is found to be the same 
as in the cr — > cr' channel in the calculation. On the other 
hand, an extra peak has been observed in the a ^ i:' 
channel at hio = 7782 eV (pre-edge peak) for ip = 0.^^^ 
This peak may come from the electric quadrupole {E2) 
transition. 

Figure 6(a) plots the peak intensity as a function of 
\Ed\ {a.t ip = —45° in the <j — > cr' channel). The intensity 
of the "main" peak increases with increasing values of 
\Ed\. It is nearly proportional to \Ed\^. 

Figure 7 shows the azimuthal angle dependence of the 
main peak intensity for G = (00|), in good agreement 
with the experiment. The same dependence as eq. (25) 
has been proposed on the basis of the symmetry of the 
scattering tensor. -'^^^ Note that the intensity ratio be- 
tween the a ^ a' channel and the a ^ n' channel is 
determined by a geometrical factor; the oscillation am- 
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co[eV] 

Fig. 5. RXS spectra for G = (00|) at V = -45° in the a ^ a' 
channel, as a function of photon energy. Top: experimental 
spectra at T = 20 K (Phase II). ^'''^ Middle: Calculated spectra 
by taking account of the direct influence of lattice distortion with 
= —0.1 eV. Bottom: Calculated spectra by taJfing account 
of the influence of quadrupole ordering 4/ states with Ef/Df = 
-0.2. 



plitude of intensity in the cr — > tt' channel is the factor 
sin^ 9 = 0.312 smaller than that in the a ^ a' channel. 
(b) G = (ftO|) with h and £ odd integers. 

The total amplitude is given by a combination of 
M(l,a;)-M(2,a;)-M(3,a;)-hM(4,w). Thus, we have 
the final form, 






i{aco)-a^)) 




M{G,uj) 

(26) 

Combining the geometrical factor in Table I to eq. (26), 
we obtain the scattering intensity. The spectral shape as 
a function of photon energy is found almost the same as 
that for the (00|) spot, so that we omit the correspond- 
ing figure. The azimuthal angle dependence is given by 



W,(G,t^)(xi|^(u;) 



C(u;) |2 (cos^ (3 cos^ V - sin^ ■^/')^ 



,(G,a;)cx-|^(a;)-C(a;)|2 

X [cos 9 sin /J(cos /3 cos tjj + sin tjj) 
— sin^(l -|- cos^ P) sin ^ cos V')]^, 



(27) 



where is determined from tan/3 = {2hc)/{£a). Figure 8 
shows the azimuthal dependence of the peak intensity for 
the (30 1) spot. A large non-resonant intensity has been 
observed in the tr — > cr' channel, and the resonant behav- 
ior is not clear in the experiment. The non-resonant 
intensity may come from the Thomson scattering due to 
the lattice distortion. We hope that the resonant behav- 
ior discussed here is observed in the a ^ n' channel, 
since the non-resonant intensity is expected to disappear 
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Lattice Distortion and Resonant X-Ray Scattering in DyB2C2 
(b)lnfluence of the 4f states 




Fig. 6. (a) Main peak intensity as a function of the crystal field 
parameter \Ej\ on the 5d states, which is calculated by taking ac- 
count of the direct influence of lattice distortion, (b) Main peak 
intensity as a function of the crystal field parameter \Ef\/ Df on 
the 4/ states, which is calculated by taking account of the in- 
fluence of quadrupole ordering 4/ states (open circles). Crosses 
represent the quadrupole rrioirients {0^2_y2) in the local coor- 
dinate frame for each sublattice. 
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Fig. 8. Azimuthal angle dependence of the RXS intensity of the 
main peak on G = (30|). 
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Fig. 7. Azimuthal angle dependence of the RXS intensity of the 
main peak on G = (00 1), in comparison with the experiment at 
T = 20 K.i^) 



in this channel. 

4-2 Influence of quadrupole ordering 4f states 

The initial state is evaluated in § 2, where the 4/ 
quadrupole moment is ordered. The intermediate state 
is evaluated by the following steps. Let the E\ transition 
take place at the origin. The complex of 4/ electrons and 
the 2p hole is assigned as the eigenstate \v) with energy 
Ev, which is calculated by diagonalizing the matrix of 
the intra-atomic Coulomb interaction. To keep the ma- 
trix size manageable, the space of 4/ states is restricted 
within the space of J = 15/2. This restriction causes 
only minor errors in the RXS spectra, since the RXS am- 
plitude contains the overlap between the 4/ states in the 
intermediate state and that in the initial state, which be- 
comes very small for the 4/ states outside the J = 15/2 
subspace. Since the photoexcited hd electron interacts 



with the 2p hole and with 4/ electrons, the complex of 
4/ electrons and the 2p hole serves as a scatterer to the 
bd electron. Thus we have the resolvent at the origin, 

( 

\hv - Hint +ir/^<i^d^.^,d^,d^, 

— Ujji^d gdi,.„^ld gldi,i\ , (28) 

where m'^ and s"^ specify 5rf states. The potential 
U^dgd^.^idgid^i includes the 5d-Af Coulomb interaction 
as well as the 5d-2p Coulomb interaction. They are ex- 
pressed in terms of the Slater integrals, which are given 
in Table 11.^^' Since the 5d-4/ Coulomb interaction is 
implicitly included in the 5rf-band energy, we eliminate 
the average of the 5d-Af interaction from the potential 
U. The energies of the 4/ states coming from the crystal 
field much smaller than other energies in the intermedi- 
ate state, and thus can be neglected. We do not include 
the crystal field on the 5d states due to the lattice distor- 
tion (the second term of eq. (15)). Within the present 
approximation, the right hand side of eq. (28) becomes a 
matrix with dimensions 640 x 640 (640 = 10 x 64), which 
are numerically inverted. 

The resolvent thus obtained is the same on four sub- 
lattices. The scattering amplitudes become different on 
different sublattices after the transition matrix elements 
from the initial state are taken into account. The scat- 
tering amplitude is found to take the same forms as eqs. 
(23). Therefore, the azimuthal-angle dependence is the 
same as eqs. (25) and (27). As regards the photon- 
energy dependence, the bottom panel in Fig. 5 shows 
the RXS spectra at ip = —45° in the a ^ a' channel 
for G = (00|). We put Ef/Df = -0.2. The spec- 
tral shape is not so different from the curve given in the 
preceding subsection, although a small hump appears 
at the low energy side. Difference is the magnitude of 
the intensity. Figure 6(b) shows the main-peak inten- 
sity as a function of |ii^/|/I?/. It increases with increas- 
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ing values oi\Ef\/Df. In the same figure, we have also 
plotted the sublattice quadrupole moment as a function 
oi \Ef\/Df. It also increases with increasing values of 
\Ef\/Df. This coincidence is plausible, since the charge 
anisotropy in the 4/ states increases with increasing the 
ordered quadrupole moment. However, its magnitude 
remains much smaller than those given by the direct in- 
fluence of lattice distortion in a wide parameter range of 
crystal fleld by comparison with the curve in Fig. 6(a). 

5. Concluding Remarks 

We have studied the mechanism of RXS at the Lm 

edge in the quadrupole ordering phase of DyB2C2. Hav- 
ing analyzed the effect of the bucking of sheets of B and 
C atoms on the 5d and 4/ states, we have constructed 
an effective model that the crystal field is acting on the 
hd and 4/ states with the principal axes different for dif- 
ferent sublattices. We have calculated the RXS spectra 
in the Ei process by treating the bd states as a band and 
the 4/ states as localized states. We have considered two 
mechanisms separately that the lattice distortion directly 
modulates the Sd band and that the charge anisotropy 
of the quadrupole ordering 4/ states modulate the hd 
band through the bd-Af intra-atomic Coulomb interac- 
tion. We have found that both mechanisms give rise to 
the RXS intensities on (00|) and (/iO|) spots with simi- 
lar photon-energy dependences and the same azimuthal 
angle dependence. Both explain well the experimental 
RXS spectra. However, it is shown that the former mech- 
anism gives rise to the intensity much larger than the lat- 
ter one for a wide parameter range of crystal field. This 
suggests that the main-peak of the RXS spectra is not a 
direct reflection of the quadrupole order but mainly con- 
trolled by the lattice distortion. To confirm this obser- 
vation more quantitatively, band structure calculations 
may be useful since the states are considerably ex- 
tended in space. This study is left in the future. 

As regards the pre-edge peak, we have estimated its 
intensity within the E2 transition. In that estimate, 
we have used the same initial state as discussed above 
and have taken account of the full multiplets of the 
/^'^-configuration for the intermediate state. The tran- 
sition matrix element has been evaluated by the atomic 
Hartree-Fock wave function. The pre-edge peak inten- 
sity thus evaluated is found to be more than three-order 
of magnitude smaller than the main-peak intensity eval- 
uated by the mechanism of the charge anisotropy of the 
quadrupole ordering 4/ states. This is inconsistent with 
the experiments, where the pre-edge peak intensity is 
the same order of magnitude to the main peak intensity. 
Clarifying this point is also left in the future. 
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